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Abstract 

A well known problem with EOP prediction is that a prediction strategy proved 
to be the best for some testing period and prediction length may not remain as such 
for other period of time. In this paper we consider possible strategies to combine 
EOP predictions made using different analysis technique to obtain a final prediction 
with the best accuracy corresponding to the least prediction error between input 
predictions. This approach can be used to improve the short-term real-time EOP 
forecast. 
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1 INTRODUCTION 



Prediction of the Earth Orientation Parameters (EOP) is a practically very important 
and theoretically very interesting task, one of the main fields of activity of operational 
EOP services. Various methods are developed to compute a highly accurate EOP forecast. 
However, a usual and well known problem is that different methods show different accuracy 
at different time intervals and prediction lengths. A method which is the best for short- 
time prediction may not be such for long-time prediction, and vice versa. On the other 
hand, a method which was proven to be the best for a testing period of time may not 
remain as such for the coming period. As an example, results of actual UT1 predictions 
are shown in Figure [H One can see that different predictions show the best results for 
different years and length of prediction. 




Figure 1: Results of actual UT1 predictions for three years. BA - Bulletin A prediction, 
AM - prediction of the same NEOS series made by the author, SI - author's prediction 
of the SLR series computed at the Institute of Applied Astronomy. Last two predictions 
were made making use of method developed by Malkin & Skurikhina, 1996 

A practical consequence is that a prediction procedure optimally adjusted for some 
time period may not remain the best for the following time period. Evidently, the main 
reason of this is that the Earth rotation is more complicated process than we are able 
describe by our forecast models. So, the question is whether we can improve our prediction 
strategy to make it more robust to unpredictable behavior of the Earth? 

In this paper we try to use combining procedures to achieve this goal. We will con- 
centrate on short-term prediction up to 5-day length as evidently the most interesting for 
practical needs, GNSS applications in the first place. 



2 TEST PREDICTIONS 

For testing purposes, nine series of everyday predictions for the period from 1 October 2006 
through 31 December 2007 making use of prediction strategy developed at the Institute of 
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Applied Astronomy (Malkin & Skurikhina, 1996). This method includes three prediction 
techniques: 

1. Least squares fitting for trend and several harmonics, 

2. Autoregression. 

3. Autoregressive Integrated Moving Average. 

Final prediction is composed of several segments computed using different methods 
and/or model parameters and merged in one continuous series using boundary conditions. 
The main model parameters that can be varied to adjust the prediction procedures are 
the following: 

1. base interval, 

2. number and period of harmonics, 

3. trend order, 

4. AR order, 

5. IMA order, 

Abovementioned nine prediction series were computed using various sets of listed pa- 
rameters aside from item 2 which remains the same for all the predictions. 



3 COMBINED PREDICTIONS 



Using prediction time series computed as described above, three combined predictions 
were computed for the period from 1 January 2007 through 31 December 2007. They are 
the following: 

CI — average of input prediction as proposed and tested by Luzum et al., 2007, and 
Schuh et al, 2008. 

C2 — best previous prediction for given length. Normally, we compute EOP prediction 
on the day of the last observed epoch ("today"). To compute n-day ahead prediction 
we examine the set of predictions made n days ago and select one that predicts today 
EOP most accurately. Then we compute our today n-day ahead prediction making 
use of the method used for computation of the best prediction made n days ago. 

C3 — best yesterday prediction. To compute this prediction we examine the set of pre- 
dictions made yesterday and use for today prediction the method which corresponds 
to the prediction made yesterday and predicts today EOP most accurately. 
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Figure 2: The RMS and maximum errors of the input and combined predictions. 



4 



Note that C2 and C3 predictions also use CI results along with the input data. The 
accuracy of the nine test predictions and three combinations are shown in Fig. [2J 

From this test one can see that all the tested combined predictions show similar accu- 
racy for the Pole coordinates, better than accuracy of input series. Marginal advantage 
of C2 and C3 can bee noticed. 

There is more complicated situation with the UT1 prediction. The average (CI) 
prediction does not provide a satisfactory result, probably because of presence of bad 
predictions (note the worst UT1 prediction in the upper left part of corresponding plots). 
The C3 (best yesterday prediction) is clearly the best for short-time UT1 prediction. 

4 CONCLUSIONS 

From our study presented in this paper we can draw the following conclusions. 

1. A procedure based on computation of multiply predictions with further combina- 
tion/choice of the best one can provide robust accurate operational EOP forecast. 

2. To achieve the best results, input predictions should be computed using the best 
methods. 
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